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Fivesquare 


In a5 x 5 array, there are exactly 50 sets of four 
squares that themselves form a square (connecting their 
centers). This agrees with the formula on page 6 of 
our issue 89 (in the article on NOSQUARE). 


The 50 sets of four squares are shown in the 
Figure on the cover. 


The numbers from 1 to 25 can be arranged in a 5 x 5 
array in (25!) ways: 


= 1551121004 3330985984000000. 


For any given arrangement, we will apply a score according 
to the following schedule: 


If all four numbers in 
a sub-square are prime 
(counting one as a prime). .. . tally 4 points 


If all four numbers in 
a sub-square are composite. .. tally 1 point 


If the four numbers in 
a sub-square have a 
factor in common .. mete tally 1 point 


If the four numbers in 
a gub-square form an 
arithmetic progression... tally 10 points 
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Test 
data 
for 
the 


FiveSquare 


problem. 


Problem: What arrangement will then produce the 
greatest score? 


= a Oe ae a 


The array represented by A(I) = I (that is, the array 
with the numbers in normal order) has a score of 34, 
consisting of 3 squares of primes, 12 squares with a factor 
in common, and 10 squares made up of composite numbers. 


Other test arrangements are given. An arrangement 
has been found with a score of 122. 


Random arrangements of the numbers in the cells of 
the array can be created by any of the standard methods 
of creating random permutations (a dozen possible schemes 
were presented in our issue number 50). A quick and 
dirty scheme is the following: 


1. Fill the array by A(K) = K. 


2. Generate a random integer in the range from 
1 to 25; call it R. 


3. Interchange element A(I) of the array with 
element A(R). 


4. Repeat steps 2 and 3 25 times, with I 
running from 1 to 25. 


In terms of BASIC, the following subroutine will implement 
the above scheme, where RND(1) returns a fraction between 


zero and one. 


2000 FOR K = 1 TO 25 

2010 A(K) = K 

2020 NEXT K 

2030 FOR K = 1 TP 25 

2040 R = INT(25*RND(1) + 1) 
2050 

2060 

2070 

2080 

2090 
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First 
Occurrence 


100000 
100001 
100011 
100026 
100039 
100077 
100117 
100139 
100429 
100529 
100777 
101587 
104815 
103698 
123998 
238521 
331809 
583895 
691845 


10000000 
10000001 
10000011 
10000026 
10000039 
10000077 
10000117 
10000139 
10000429 
10000529 
10000777 
10001117 
10003669 
10004929 


10008539 
10014765 
10014283 
10027953 
10017790 
10014944 
10128033 
120184 33 
15418423 
31789816 
41279617 
70319091 
83099616 
88169426 


Distrib- 
ution 


901 
8031 
30179 
68826 
113017 
148529 
158301 
139609 
103073 
64549 
35819 
17542 
T7451 
2879 
986 
271 
102 

a4 

2 


9001 
117092 
580092 

1721776 
3777752 
6659395 
9774369 
12210996 
13165051 
12418609 
1037 3834 
7704 328 
5143332 
3107025 


1705741 
859351 


Problem 280 (The Six-Digit Algorithm) in issue number 92 called for 


partitioning a 6-digit number into two 3 


together to produce a new 6 


-digit parts; multiplying those parts 


-digit number; and continuing this process until 


it degenerates to zero, and counting the number of stages. 


Harry L. Nelson, Livermore, California, calculated all the results for 


both 6- and 8-digit numbers, as shown here. 


To: Popular Computing: 


Popular Computing number 92, page 3, refers toa 
Professor Donald Knuth who wrote a book called The Art 


of Computer Programming, Vol. 3, that “omits bubble 
completely." 


I wonder which Professor Donald Knuth and which 
book you are referring to, since my copy of a book with 
the same name has the following entry in its index: 


Bubble sort, 106-111, 134-135, 141, 
223-224, 244-246, 352-353, 360, 
379, 386, 388, 644, 


Of course, something that is mentioned on only 22 pages 
is easy to miss. 


The method described in that book also bears a 
curious resemblance to the flowchart given on page 7 of 
PC92, except that the book method is slightly more 
efficient and the program is shorter. 


Well, I never did think many people could under- 
stand the inscrutable language of The Art of Computer 
Programming, but I wish there were a way to spread the 


word that it does have an index. This index can be 
used (as Jonathan Swift said) by people who wish to 
pretend they have read the book. 


Cordially, 
Donald E. Knuth 


Actually, our blunder in misquoting the master is, 
indirectly, a compliment. We don't have a copy of 
Knuth, Vol. 3. Editors can obtain almost any book 
in their field for the asking; publishers of books 

_ could never risk refusing a request for a review copy. 
There is one big exception; namely, Addison-Wesley and 
the three volumes of Knuth--if you want them, you pay 
full price, no matter who you are, We will probably 
buy Volume 3 when Volume 4 comes out. 


The blunder was unintentional, but its effect-- 
a nice letter from Prof. Knuth--was delightful. 
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On the facing page there 1s empirical data of 
16 observed values (points 7 through 22). You are to 


predict, 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


as well as you can, the time for point 29. 


You will, of course, want to plot the data. 

Will the nature of the problem be altered if 

you change the scale of the dependent variable? 

Will you try a least squares curve fit? if 

so, what degree seems appropriate? 

Could you use the method of undetermined I 
coefficients? If so, what degree of poly- 

nomial would be appropriate? 

Could the problem be solved graphically, simply 

by making a large and accurate plot of the data? 

Just how accurate is the given data? What 

confidence do you have in the extrapolated » 
result? In other words, how many digits of 

your result would you be willing to guarantee? 

If you use the method of least squares, then you 

should agree that there should be just one correct 
result and everyone in your class should achieve 

the same result (to the same precision level). 


Is this statement acceptable? 
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10:27:20 


10:27:48 
9 10: 28: 37 
10 10:29:50 
The times shown in the 
11 10: 31:28 
second column are of 
12 10: 33:53 
succeasive observations, 
13 10: 37:00 
in hours, minutes, and 
14 10:41:15 
seconds, of a computing 
15 10:46: 30 
experiment conducted 
16 10:53:18 
one morning. We seek 
17 11:01:15 
an estimate of the time 
18 alabp ibe aS) 
for the 29th observation. 
19 11:23:00 
20 11: 37:00 
el 11:53:05 
22 12:12:40 
Lg 
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CURVE FITTING and EXTRAPOLATION 


al 


Progressive Mantissas 


Consider the following algorithm. Start with any 
number, say 2, and call it A. Take its square root 
and add the mantissa of the square root to A. 


Repeat this process indefinitely. Keep a count, N, 
of the number of terms and calculate the ratio N/A = Q, 
Thus, we develop a sequence that begins: 


N Progressive sum 


1 241421356 
2 2.96798754 
3 3 .69077236 
4 
5 


4 61191066 
475944661 1.0505 


Now, an average mantissa should be around .5, so 
the values of A should increase half as fast as N, and the 
ratio Q@ should tend toward 2.00. But it doesn't--why? 


For various starting values, we have the following 
progressive sums and ratios at 1000 terms: 


Starting A Progressive sum 


«345590 
.561512 - 
- 559193 
.210541 
.054038 
. 748623 
. 346720 
. 371066 
219195 
.6457 32 
556785 
-984940 
.925925 
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e Problenn Solution 


Problem 195, REVERSE, in issue 55, was this: 


Start with all the natural numbers, and 
reverse every group of K = 2, producing the 
sequence: 


2, i, 4, 3, 6, 5s 8, Ts 10, 9, nos ii, oa. 


Print the first set of 2 and delete them from 
the sequence. Increase K to 3. Reverse every 
group of 3, to produce: 


Bese is Oso, ae, 9p 105° 135 1a ee 


Print the first set of 3 and delete them. 
Increase K to 4, and reverse each group of 4, to 
produce the sequence: 


@ 12, 5, 8, 7, 14, 13, 10, 9, 16,... 


Increase k, and so on. Calculate the first 
1000 numbers to be extracted by this process. 


The problem lends itself nicely to calculation in BASIC. 
The logical steps are these: 


i. Create an array, B, of dimension 2000. 
(and an auxiliary array, C, also of 2000 terms) 


2. Fill array B with consecutive numbers; that is, 
set B(I) = I. 


Set K = 2, 
Reverse all groups of K terms. 


Print the first group of K terms; namely, 
Bil), Biv. B(3), Nev BK) 


Shift the array left K terms. 
Increment K by l. 
Go back to step 4 and repeat. 


In terms of Applesoft BASIC, some of these steps can be 
carried out as follows: 


FOR I = 1 TS 2000 
C(I) = B(T) 
NEXT I 


FOR J = 1 TH Te STEP K 
FOR L=1 MK 

B(J+K-L) = C(J+L-1) 

NEXT L 

NEXT J 


RETURN 
FOR I= 1 TOK 
PRINT I, B(I) 
NEXT I 

RETURN 

FOR I = 1 TW 1900 
B(I) = B(I+K) 
NEXT I 


RETURN 


The first four successive groups that are printed are: 
1 2 16 
2 1 
3 
4 
5 


and the next seven groups are these: 
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Db F O 


Examination of these, and successive groups, shows a 
pattern. If we list, for each group, the number of the 
largest term and the number of the smallest term, we get 
the accompanying table. 


Will this pattern continue indefinitely? That is, 
can we say with assurance that when we reverse groups of 
1002 terms, the 1002 numbers that are extracted and printed 
Will have their largest number at the 100lst term and 
the smallest number at the 502nd term? 


Incidentally, the group for K = 45 will print out as 
follows: 


1352 1205 1218 1149 1150 1143 1144 
1108 1055 1062 1045 1046 931 932 
896 889 890 821 834 795 828 
832 857 874 909 g4e 961 968 
978 981 1070 pip tal “plik 1115 1138 
1210 1229 1236 1245 1282 1353 


Thus, the 1000th number to be printed is 1062. 


2 1 2 26 25 14 
3 1 2 27 1 14 
4 1 2 28 1 14 
5 1 2 29 1 14 
6 5 4 30 29 16 
" 1 4 31 1 16 
8 1 4 32 1 16 
9 1 4 33 1 16 
10 9g 6 34 33 18 
11 1 6 35 1 18 
12 1 6 36 al, 18 
13 3 6 37 1 18 
14 13 8 38 37 20 
15 1 8 i 39 il 20 i 
16 1 8 | 40 1 20 
17 1 8 j 41 1 20 
18 17 10 ve) yy 22 
19 1 10 43 1 22 
20 1 10 44 i, 22 
el 1 10 45 1 22 
22 21 12 46 45 24 
23 1 12 47 1 24 
24 1 12 48 1 24 
25) 1 12 49 1 24 
Me / y 50 4g 26 
This is K eA | 
\ 
ie group K, this is the number of the \ 
term with the highest value 


In group K, this 1s the number 
of the term with the lowest value. 


A pattern observed in the REVERSE problem. 
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Computer Economics 


The following square roots of integers: 
W 2.64 prelees 
31 ead) 6 4 3 
37 8 
42 8 
58 : 1 
76 1 
78 3 
D5 9.74 
have a feature in common; namely, for each of them, the 
4th decimal place is a 7. We might expect that, on the 


average, one number in ten (not necessarily integers) would 
have a square root that exhibited that characteristic. 
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2 
O 
3 
7 
O 
4 


6 
4 
i 
9 
6 
9 


7 3 


But restricting our attention to the square roots of 
integers: suppose we want a long list of such numbers. if 
the long list is only a few dozen, we would consult a printed 
table of square roots (although such things are getting 
scarce); put a straightedge along the 4th decimal column; 
and read off results by eye. Including the time it takes 
to record the numbers we seek, this method should proceed 
at the rate of about one per minute. 


On the other hand, if we need a million such numbers, 
we would have no hesitation in writing a computer program 
to sift them out. 


Where is the dividing line between those two extremes? 
There are many tasks that could be performed manually (or 
semi-manually) by clerks, or that could be automated by 
computer; how does one make a rational decision between the 
two courses of action? The cost of the manual approach 
is usually weighed against the cost of writing, debugging, 
and testing a computer program, plus the cost of CPU time 
(although, as in the case of this problem, the cost of the 
computer itself may be negligible). 
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Suppose we make the following assumptions: 


(1) Programmers cost four times as much as 
clerks. 


(2) The problem can be analyzed, coded, 
debugged, tested, and production begun, 
in four hours. 


If we go the computer route, then after the program is in 
production, the cost per unit result will approach zero 
the longer we run the program. If we go the clerk route, 
the total cost will be less up to 16 clerk-hours, so 
whatever amount can be done in that time is the cutoff 
value. The trivial task we have outlined would best be 
done by two clerks working together; in 16 clerk-hours we 
might expect to get 500 correct results. 


The problem involved here suggests two further 
problems: 


(1) Considering the square roots of successive 
integers, what is the longest string of 
consecutive integers that one can find for 
which each square root does NOT have 7 in 
the 4th decimal place? 
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(2) For any numbers, it is to be expected that 
about one in ten square roots will have a 
7 in the 4th decimal place. The 
accompanying flowchart suggests a way to 
demonstrate this. The functions on the 
flowchart are those of BASIC; namely, 
INT for the greatest integer in a number, 
and RND(1) for a new random number between 


zero and one. [eal 


Cycle Lengths 


It is well known that the low-order digits of powers 
of 2 cycle as follows: 


O46" 5.76 


Nw by 4's 
| 
\W_-—— by 20's 
a, eer 
by des 


| | 
2 -by 500's 


re 


| | 
Mee Pe be by ee cer 


A similar pattern is exhibited by powers of 3: 


the unit's digit cycles by 4's 

the low-order two digits cycle by 20's 
the low-order three digits cycle by 100's 
the low-order four digits cycle by 500's 
the low-order five digits cycle by 2500's 


5 oatel EXe) Cyaic o 5 


Powers of 5 have a distinctly different pattern. 
After the first power, all powers end in 25, and the cycle 
length for the R low-order digits is as follows: 


The table shown here 
gives the first 25 
owers of 2 and 3 
that is, their low- 
order three digits) 
and sufficient powers 
of 5 to demonstrate 
the repetition in the 
low-order digits. 
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-..one might conjecture that the cycle length for the seven 
low-order digits of powers of 5 is 32. But conjectures 
like that are dangerous. 


And what of powers of other numbers? Powers of 7 
seem to have the same pattern in their cycle lengths as 
powers of 2, but powers of 11 are quite different: 
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It would seem that there is room here for a modest 
piece of research, which could be done on any machine in 
any Sea oe (the results given above were obtained in 
BASIC). 


In the book Programming in BASIC for Personal Computers 
(David Heiserman, Prentice-Hall) there appears the following 
amazing paragraph: 

The function EXP(n) simply raises any number n 
to the eth power, where e is a constant value very 
close to 2.71828. The same sort of function can 
be carried out by means of the BASIC formula 
n#2.71828. But, of course, EXP(n) is far easier 
to use and it runs much faster on the computer. 
Again, this is one of those functions normally 


reserved for engineering and scientific applications. 


Is there any value of n for which this weird distortion 
could be true? That ig, for what value of n, if any, 
is it true that: 


Doubtless this question can be settled analytically, 
but 1t 1s posed here as a computing problem. 


On the assumption that there is a value of n to be 
found, set up a search procedure (bracketing would be 
ideal) for it. If your search fails, can you then 
condlude that there is no value of n that satisfies? 
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